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ABSTRACT 


This  report,  prepared  in  two  parts,  deals  with  products  and  quotients 
of  random  variables.  In  Part  I,  the  distributions  of  quotients  of  indep^  dent 
random  variables  are  considered.  In  Part  II,  the  distribution  of  the  product 
of  two  (not  necessarily  independent)  normally  distributed  random  variates  is 
investigated.  The  tables  of  this  distribution  are  given  in  the  Appendix. 
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I.  ON  Tim  QUOTIENT  OF  HAWPQH  VARIABLES 


1.  Introduction 

In  a  study  of  distributions  of  products  and  quotients  of  random 
variables  it  is  sometimes  necessary  to  determine  possible  component 
distribution  when  the  composite  distribution  is  known.  Formally,  at 
least,  this  involves  a  study  of  linear  integral  equations  of  the  first 
kind.  For  the  quotient,  in  particular,  suppose  x^  and  x 2  are  independent 
random  variables  with  f(x^)  and  g(x their  respective  density  functions. 
Setting  y^  -  x\/x2  ^2  "  x2*  density  function  for  the  quotient 

oo 

has  the  fora  y?(y^)  -  /  f(y^y2)  g(y2)  J  y2  |  dy2»  A  derivation  of  this 

-oo 

formula  as  well  as  a  general  discussion  of  results  in  this  area  is  given 
in  £  20 J  „ 

A  number  of  authors  22  f  2  3  f  51  f  32f  33,  3\  and  others  j 
have  studied  this  problem  in  the  case  where  the  variates  are  assumed  to 
be  identically  distributed.  Their  techniques,  which  can  be  called  more 
or  less  "classical",  involved  the  use  of  various  transform  theories  — 
Mellin,  Fourier,  and  others  —  but  little  was  done  to  develop  a  general 
theory  for  the  above  equation.  It  was  the  original  intent  of  this  paper 
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lanuscript  re  loused  by  the  authors  Dcccnltr  19^4  for  publication 
as  an  ARL  technical  report . 


to  view  this  equation  in  operator  form,  i.e. ,  F  g  (assuming  f  given), 

where  F  is  a  compact  operator  on  an  space  and  attempt  to  develop  a 
theory  for  linear  integral  equations  of  the  first  kind  from  the  linea” 
operator  point  of  view  that  would  be  applicable  to  probability  density 
functions.  We  have  not  been  too  successful  as  far  as  fruitful  results 
are  concerned.  Some  of  the  difficulties  are  discussed  in  Sec.  3  along 
with  a  possible  application  of  a  theory  that  is  presently  being  developed 
for  pairs  of  operators  acting  between  Banach  spaces. 

In  Sec.  2  we  generalize  a  few  known  results  assuming  identical 
distributions  as  mentioned  above  (still  from  the  classical  point  of  view) 
and  discuss  an  existence  theorem.  Essentially,  the  procedures  and  tech¬ 
niques  of  Laha  have  been  followed  here.  In  Appendix  A  two  theorems  re¬ 
sulting  from  a  side  investigation  are  presented  -  one  is  apparently  not 
new  but  the  proof  seems  particularly  simple. 
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2.  The  Quotient  of  Independent  Random  Variables 


We  give  here  a  few  results  and  generalizations  of  known  results 
involving  the  quotient  of  independent  random  variables.  In  parts  of 

*t  ii 

Seo.  2.1  and  Seo.  2.2  identically  distributed  is  also  assumed. 

2.1.  A  Necessary  Condition  for  a  Solution  to  Exist 

Suppose  x^  and  x^  are  identically  distributed  random 
variables  over  the  real  line  with  density  function  f.  Under  the  trans¬ 
formation  y^«  X\IX2  an<*  ^2*  X2  •  8ee  fco]  »  the  density  function  for  the 
quotient  x^/xg  *8  <?iven  by 


(A) 


■  / 

-oc 


f(y2)  |  y2|  dy2* 


Hote  then  that  ■f  f(y2/y1)  f(y2)  |y2  |  dy2-  For  y^  o  ,  letting 

'  ^  -00 

y^/y^-  w  so  that  dy^*  y^dw  this  beoomes 


r(w)  f (yiw)  I yi  1 1 w  1 1 yi  I 

-00 


dw 


00 

-  |y i\2  J"  f(w)  |w  |  dw 

-00 


-  y 
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Note  also  that  ^(0)  =  f(O)  E  |"  J  y^  jj  .  Viewing  j  as  given  in  (A)  above 
we  h  ive  proved  the  following. 

Theorem  1»  Under  the  conditions  above,  if  (a)  has  a  solution  then  for 

0  it  is  necessary  that^^^-^  =  y^^(y^);  if  not,  no  solution  f 

exists.  Moreover,  the  "search"  for  possible  solutions  can  be  narrowed 

down  to  those  density  functions  f(x)  such  that  ^0)  °  f(o)  E  j"  |  x  |J  . 

2 

It  is  interesting  to  note  in  this  result  that  y^  ^(y ^ )  is  actually 
the  density  function  for  the  random  variable  l/y^.  Note  further  that  (a) 
cannot  be  solved  if  it  is  assumed  that  (y^)  is  the  normal  density  function. 
Referring  to  (A)  again,  suppose  a  symmetric  solution  is  desired,  i.e., 
f(z)  =  f(-z).  The  equation  then  has  the  form 


00 


oo 

f(y1y2)  f(y2)  y2  <*y?. 


It  is  clear,  first  of  all,  that  f  symmetric  implies^  must  also  be 
symmetric.  Also  if  (j'  is  given  symmetric  then 


OD  00 

J  I  y2  |  f(-yxy2)  f(y2)  dy2  =  J  \  *2 1  f!^i>r2^  f^2  iy2 

-OD  -00 

which  implies  that  =■  fCy^y^  a.e.  Honce,  in  the  a.e.  sense 

at  least,  we  have  the  following  for  (a)  above. 


Lem*2__2t  is  symmetric  if  an  only  if  f  is  symmetric. 
Changing  the  form  of  (B)  above, 

oo 

f^yl-  “  2 f  y2*  f(y2^  y2^  f^yly2^  dy2  * 

o 

Letting  y^2  f(y?)  3  this  bec0Bua 


oo 

2  f  g(y2)  y2  f(y1y2)  dy2 

O 


BO 


OO 

that  for  yx  >  0,  y^  =  2  f  g^y2^  yl^  y2^  f(yiy2^  ^2 


*  i»e. , 


oo 

e(yi)  -  y^fCy^  -  ?  J  g(y2)  g(yxy2)  <*y2 


This  is  a  form  studied  in  21j  .  Solving  the  above  equation  for  g  will 
also  furnish  a  symmetric  solution  to  the  original  equation. 

It  should  be  mentioned  that  Fox  |^2lJ  carries  out  an  analysis  on  the 
above  form  for  L  'otoo)  functions  using  Mellin  and  Fourier  transform  theory. 
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2.2.  Some  Generalizations  Using  Both  Diatrlbution  and  Density 
Functions 

Laha  considers,  in  particular,  an  integral  equation  of 

the  form  (A)  above,  where  y  ,  the  quotient  of  two  independent,  identically 
distributed  random  variables,  is  assumed  to  follow  the  Cauchy  law.  The 
general  technique  is  to  use  distribution  functions  and  Fourier  transforms, 
the  distribution  functions  assumed  to  be  everywhere  continuous  to  the  right. 
This  is  a  more  general  approach  since  the  distribution  function  alwayB 
exists. 

The  distribution  function  F  for  the  random  variable  y  is  said  to  be 
symmetric  (about  0 )  in  case  F(y)  =»  1  -  F(-y-O). 


Lemma  li  Given  the  ruuidom  variable  x  with  distribution  function  F(x) 
symmetric  (about  0) ,  the  distribution  function  G  of  jx  |  is  given  by 


2  P(x)  -  1, 

0, 


x  >  0 

el sewhere  . 


1  roof :  for  ft  >  J, 


ha; 


Pr 

- 1 

X 

A 

23  a 

*  Pr 

-a  <  x  <  a 

L  ~  J 

F(a)  -  F(-a-O)  -  2F(a)-l 


since  F  is  symmetric  and  the  result  follows. 

The  importance  ol  the  lemma  is  that  the  distribution  function  of  x  can 
be  determined  knowing  only  that  of  I  x  |  and  we  shall  be  able  to  relate  this 


to  the  distribution  of  In 
following. 


.  Before  doing  this  however  we  prove  the 
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Lemma  2t  Let  u  be  a  random  variable  following  the  Cauchy  law,  i.e., 


f(u)  -  - -  --  .  Then  z 

>7(l«i2) 


cot  u  has  a  uniform  distribution. 


Proof i 


so  that  the  density  function  for  z  is  equal  to  ,  0  <  z  <  77"  and  zero 
elsewhere.  Using  the  above  lemma,  since  z  has  n  uniform  distribution,  and 
since  we  knew  its  closed  form  characteristic  function*  see  for  example, 

and  moreover  ainoe  the  characteristic  function  of  a  function  of  a  random 
variable,  g( w)  say,  is  the  mean  value  of  e*  *  ^  ^  we  can  evaluate  the 
following  integral, 


lt/r/2  aln  urjz  ■ 
*  \rr/2 


Theorem  3 i  Let  x  and  y  be  independent,  identically  distribut’d  random 
variables.  Let  z  -  x/y  and  G(>),  the  distribution  function  of  z,  be 
symmetric  about  0.  Suppose  further  that  the  square  root  of  the  charac¬ 
teristic  function  of  In | z  |  is  absolutely  integrable.  Then  F(x),  the 
distribution  function  of  x  (and  y)  is  absolutely  continuous  and  has  a 
continuous  density  function  f(x)  «  P'(x)>0. 
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Proofi  By  a  result  of  Lana  22J  ,  F(x)  is  synraetric  about  the  origin. 

[2F(x)-l ,  if  x  >  0 


Consequently  the  distribution  function  for  |  x  |  is  G(  x  ) 

by  the  above  lemma.  Let  *  in|z  |j  =  (f ^ (t).  Then  3ince  In 

In  1  y  I  we  have 


0  ,  otherwise 
z  I  »  In  I  x  I  - 


/«( *> '  ■  r.(t) 

andh.no.  |^(t)|  . 


By  assumption 


ao 

■/ 

-OD 


<  00  so  that  the  '  aaracteristic  function 


of  ln|z  J  ,  <^(t) ,  is  absolutely  integrable.  By  a  theorem  of  Loeve  |  24j 
the  distribution  function  of  lnjx  |  is  absolutely  continuous  and  has  a 
continuo  s  density  function,  but  since  F(ln|  x  ( )  =  G(|x|)  it  follows  that 
|x  j  has  an  absolutely  continuous  distribution  function  and  a  continuous 
density  function.  From  above,  then,  so  does  x. 

It  is  known  that  if  z  follows  the  Cauchy  law,  the  characteristic 
function  for  ln|z  |  is  sech(/7't/2  )*  a  function  which  has  a  finite  integral 
over  the  real  line  so  that  the  above  result  holds  for  this  particular 
distribution. 
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Theorem  4»  To  the  assumptions  in  theorem  3  above  add  that  arcot  z  has 


a  uniform  distribution.  Then  f,  the  density  function  for  x,  satisfies 

the  integral  equation 
oo 

f(y)  f (wy)y  dy  -  - r  ,  where  k  is  a  constant. 

1+w 
o 


Proof i  Our  assumptions  imply  that  f  is  also  symmetric  about  zero  and 
hence  wo  recognize  the  above  integral  as  —  g(w),  where  g  is  the  density 
function  for  w  -  x/y.  Let  u  -  arcot  w.  The  density  function  for  u, 
h(u)  say,  is 


b.(u) 


““  *  0  <  u  <  rr 

0*  elsewhere, 


so  that  for  the  distribution  function  H(u)  we  have 


cot  a 


H(u)  =■  Pr|^u<aJ-  1  -  Pr  ^  w  <  cot  aj  -  1  -  j  g(w)  d(w) 

.00 


%' 0  ^ a  s r{’ 


Hence 


oot  a 


/ 

-00 


\  g(v)  dw  -  ~  ~  {rf-  a) 


cot  a 


cot  a 


But 


J  dw  =  (rf-z)  80  that f  — ~  dw  =  {ff -  a), 

-OD  -00 
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1  k 

Thus  —  g(w)  - r  and  the  result  follows.  This  can  be  generalized 

2  1+w 

as  follows. 

Theorem  Given  the  hypothesis  of  Theorem  5»  let  w  =  x/y  and  h(u),  the 
density  function  for  u  -  arcot  w,  vanish  outside  the  interval  j^O , 77" J .  Then 
the  density  function  f  for  x  satisfies  the  integral  equation 


oo 

/ 


f (y)  f(wy)y  dy 


k  h(cot',^w) 
1+w2 


,  k  a  constant. 


Proof i  Following  the  model  and  notation  of  the  laBt  proof  we  have 


co 

/. 


cot  a 


\  g(w)dw  - 


1  -  H(a) 


,  o  <  a  <  rr. 


f  2  1 

We  need  a  function  !i  such  that  £(oot  a ''esc  a  -  —  h(a), 


i.e.  i(w)  =  \  h^COt~  *1  m  Consider  \ 

1+w 


cot  a 


h(oot  Xw) 


/‘U.Wr.Y-S - ZJL  Ay 

'  2 


-00 


1+w 


Letting  y  *=  cot  ^w  this  becomes 


_  lj  h(y)dy  .  -J  h(y)dy  .  -  |^H(/T  )  -  H(a)j  -  j^l  -  H(a)j 

2Jrr  a 


and  the  result  follows. 
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}.  Some  Coaments  on  the  General  Problem 


It  has  already  been  noted  that  the  density  function  for  the  quotient 
of  two  random  variables  x^  and  x 2  with  density  functions  f(x1)  and  g(x2) 
respectively  has  the  general  form 

oo 

^(y^  » J  |  y2  |  f(yxy2)  «(y2)dy2 

-00 

where  y^  and  y2«  x2.  Writing  the  kernel  k(yL,y2)  -  |j2|  fCy^) 

this  has  the  form  of  a  linear  integral  equation  of  the  first  kind  and  in 
operator  notation  can  le  formally  written  as  Kg  -  ^ .  It  appears  somewhat 
difficult  to  determine  the  proper  domain  and  range  for  this  operator  so 
ae  to  apply  directly  to  probability  functions.  The  set  of  density  functions 
in  C  £a,bj  or  1?  [  a,bj  for  example,  does  not  form  a  linear  space.  If  the 
equations  could  be  modified  to  consider  distribution  functions,  addition 
can  be  defined  asF  +  G*F«G  where  *  means  convolution  but  scalar 
"multiplication"  appears  to  move  one  off  the  intersection  of  the  unit  ball 
with  the  positive  cone  in  the  Banach  space  under  consideration. 

Apparently  then, the  analysis  should  be  done  on  some  other  space 
(as  far  as  solutions  are  ooncemed)  and  a  second  analysis  done  to  determine 

whether  the  solution  or  which  of  the  solutions  ere  Hensitv  functions.  If, 
for  example,  the  kernel  k  vanishes  outside  the  square  [a,bj  x  [a,b]  »  £ 

2  2 
and  k£  L  (E),  then  the  operator  A  acting  on  L  is  compact  (with  range  in 

2 

L  ).  The  theory  of  compaot  operators  could  possibly  be  extended  so  as  to 
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apply  to  aquations  of  the  type  needed  here.  To  this  end  it  should  be 
mentioned  that  S.  Bimbaua,  at  the  University  of  Colorado  and  Martin* 
Denver,  is  presently  developing  a  theory  for  pairs  of  operators  acting 
between  Banach  spaoes.  This  theory  appears  to  have  sons  applications  in 
this  area.  We  give  a  brief  discussion  here  as  to  the  type  of  results  to 
expect.  They  will  be  called  pretheoreas. 

We  will  be  considering  the  spaoes  L  •  Lp  [o,l]  and  C  -  C  [o,l^  , 
1  <p  Coo ,  and  an  integral  equation  of  the  above  fora  Kf  •  g,  KsL — * C. 
With  proper  restrictions  on  the  kernel  k(s,t)  determining  K,  K  will  be  a 
continuous  operator  I  vis.  k(s,t)  continuous  in  s  for  every  t  and 
J *  j  k(s,t)|P  dt  C oo  for  every  s  £  £o,lJ  »  wher*  p  ^  p»  •  1  ]•  I*t 


R  t  C— *LP  be  an  imbedding,  l.e. ,  for  y  £“  C,  Hy  .  £yj  £“  LP,  ([,  J  denotes 
an  equivalence  class).  Then  K  is  continuous  and  1-1  and  hence  R~*  -  S 
exists  as  a  closed  operator  S  t  Lp- — >C  and  the  domain  of  S  is  the  set 
^£y]  |  ^  y  £"  [y]»  y  c  cj  .  Using  the  above  theory  for  the  pair  (S,K) 


the  following  can  then  be  proved. 


Pretheoreni  If  there  is  A  (complex)  suoh  that  (S  K)  is  1-1,  the 
range  of  (S  -  A  K)  is  C,  and  such  that  ||  S(S-  A  *)”*  ||  <  1  then  K 
restricted  to  the  domain  of  S  is  1-1,  has  range  C,  and  has  a  continuous 
inverse  defined  everywhere  on  C. 

It  should  be  noted  that  SMQK other  things  this  is  a  uniqueness  theorem. 
Applied  to  the  problem  considered  in  this  paper  the  solution  would  have 
to  be  examined  to  determine  whether  or  not  it  is  a  density  function.  Prom 
the  form  of  the  solution,  however,  it  appears  that  this  may  be  a  difficult 
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problem  but  it  has  not  yet  been  investigated.  It  also  appears  promising, 
using  the  new  theory,  that  it  will  be  possible  to  characterize  the  null 
space  of  K  so  that  something  can  be  said  when  multiple  solutions  are 
involved.  Although  the  above  result  was  stated  for  the  unit  interval, 
it  can  be  extended  to  more  general  settings. 
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II.  NOTES  ON  THE  PRODUCT  OF  TWO  NORMALLY 


DISTRIBUTED  RANDOM  VARIABLES  . 


1.  Introduction 


L«t  X.  and  X  follow  a  normal  bivariate  probability  density  function, 
p.d.f. ,  with  expected  values^,  fx^%  standard  deviations,  O^,  <J^,  and 
coefficient  of  correlation,  |0.  Several  forms  of  randcm  variable  products 
■ay  be  considered;  two  of  which  are  the  normalized  product  Z  -  (X^-  /<^)  * 
/*2 )/ °1  °2  811(1  Pro<*uc*  ^  -  X^X^/o^  6~2>  The  latter  presents  far 
greater  application  in  that  families  of  noimal  random  variables  X^/d^  and 
X^/d^  may  be  characterized  by  the  statistics  v^  « i  •  1,  2*  These, 
of  course,  are  the  reciprocals  of  the  respective  coefficients  of  variation. 

The  joint  p.d.f.  of  the  normal  random  variables  X^/cr^  and  X?/o^  iB 

1 


exp  <  - 


(¥i) 


2(1  -p2) 


27 r 


V 


1  - 


With  the  transformation  W  «  X^/d^,  Z  -  X  <5^,  the  marginal 

p.d.f.  of  Z  may  be  derived  from 


<fK  *) 


/ 

-oo 


1 

exp  j 

(  1 

(w  -  vx)2  - 

■  2p  (w-^Xj-Vg)  +  (J-  v2)2 

! 

f  2d-(02) 

27T  V 

1  1  -  /°2  M 

dw 


(ii-D 
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The  p.d.f.  (H-1)  may  be  expressed  as  ^(z)  -  IjU)  -  I2(z)  where 


and  12(2)  is  the  same  function  defined  on  (-00 ,0).  After  the  substitution 
w  =  -w  into  the  marginal  p.d.f.  ^(z)  may  be  expressed  as 


2(1  -(02) 


(w+v^^)2-  2|0  (w+v1)(-  +  v2)  +  (fj  +  v2)2  | 


dw  ; 

(11-2) 


and  by  expanding  these  exponents  and  regrouping  terms,  (II-2)  becomes 
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Several  speoial  cases  may  now  be  examined, 
p.d.f.  of  the  "normalized"  product  is  obtained, 
case  : s  < 


When  v2  “  Y2  “  0*  *he 
The  p.d.f.  of  Z  in  this 


(II-4) 


where  is  a  modified  Bessel  function  of  the  second  kind  of  zero 

order  possessing  a  singularity  at  z  -  0.  The  product  of  two  independent 
"normalized"  variables  by  (II-4)  reduces  to 


f(i)-^«0W  ,  (II-5) 

a  result  shown  in  [  5  ]  and  [  6  ]  . 

The  non-central  product  Z  -  X^X^/cr^  in  which  each  variable 
is  characterized  by  its  respective  reciprocal  of  the  ooeffioient  of  variation, 
v  4-  0,  has  undergone  extensive  study.  As  yet,  however,  no  satisfactory 
method  of  obtaining  numerical  results  for  the  cumulative  distribution 
function  of  Z  has  been  derived  for  all  parameters  values  of  p  ,  v^,  and  v^. 

The  analysis  by  C.  C.  Craig  [  5  ]» [  6  ]»  ia  perhaps  the  most  notable  con¬ 
cerning  this  product. 


J.  Wishaxt  and  M.S.  Bartletts  The  Distribution  of  Second  Order  Momenta 
Statistics  in  a  Normal  System;  Proceedings  of  the  Cambridge  Philosophical 
Society,  Vol.  28,  19?2,  PP»  ^55-^59* 
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In  an  effort  to  simplify  any  numerical  calculation,  Craig  reformulated 
( II- 3 )  as  an  infinite  series.  The  ooah  function  in  ( 1 1—3 )  Bay  bt  expanded 
so  that  it  is  possible  to  write  (f{  i)  -  I^(z)  -  I^  ( z )  where 


(II-6) 


and  I2(z)  is  the  integral  of  the  same  function  over  the  interval  (-00,0). 
The  infinite  series  expression  is  derived  by  substituting 


y  into  I  (*)  and  I2(z). 


Under  this  transformation,  I^(z)  becomes 


u 
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aeries  in  powers  of  u  for  all  u,  u^  0.  Thi6  expression  is  simplified  to 
some  extent  by  substituting 


r-1  1  V' 

In  the  expansion,  the  coefficient  of  u  ,  r  >  1,  is  — f  /  (R^R^y) 

_  r 

in  which  ^  ^  (•),  the  confluent  hypergeometric  function  of  order  r,  is2 

r 

v  ,  ,  8iV  tHiV>2 

/  (R-.R0y)  ■  l  +  i  +  7o\  ■v  / A  +  •  •  •  » 

1  2  r+1  (r+2)  2  2!  (r+})U,5! 


with  (r+k)^  =  (r+k)(r+k-l)  .  .  .  (r+1). 

By  this  expansion  and  a  similar  expansion  for  I2(y),  the  p.d.f.  of 

p 

Y  o  XjX2/o^  ff*2  (l  -  p  )  may  be  expanded  in  an  infinite  series  involving 
confluent  hypergeometric  functions  and  powers  y,  v^,  and  This  series  is 


<pw 


€ 


exp 


-1 


2(l-(0  ) 


v/+  v22-  2?[y  +  v^] 


I]  (W)Ko(y) 


a- 

♦  («i2+  »22)  ^H1H2y)K1(y)  +  (ri4+  R24)  (RiR2y)Vy^ 


(R1R2y)K}(y)  +  .  .  . 


(n-e) 


wherei  K^(y)  =  the  Bessel  function  of  the  second  kind  of  the  i^  order 


These  functions  are  disqussed  in  detail  in  Whitaker,  E.T,,  and  P.N. 
Watson,  A  Course  in  Modem  Analysis.  Cambridge  University  Press, 
Cambridge,  1958* 
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and 


(RjHjj) 


in  which  I  (-  )  is  the 
J 


Bessel  functijn  of  the  first  kind  of  the  j-  order. 

When  p  a  0,  the  p.d.f.  of  2  .  X^X^/cr^  0^  is  the  simple  Bessel 

function  expressed  by  ( II— 5 ) • 

Craig's  result  has  unfortunately  proved  to  be  of  little  use  compu¬ 
tationally  for  it  may  be  shown  that  for  large  and  (a  frequent 
occurrence  in  engineering  studies)  the  series  expansion  converges  very 
slowly;  in  fact,  for  v^  and  v?  as  small  as  2,  the  expansion  is  unwieldy. 

L.A.  Aroian [  1  ],[  2  J  took  up  the  problem  of  convergence  in  Craig' 3 

series  expansion.  Using  Craig's  notation,  he  showed  that  as  v^  and  v^ - ► 

the  p.d.f.  of  Z.  XiX2/^<Jl  °2  aPProache8  the  normal  p.d.f.  In  addition, 
he  demonstrated  that  the  Type  Hi  function  and  the  Cram-Charlier  type  A 
series  afford  excellent  approximations  to  the  distribution  of  Z  whence  =  0. 

The  characteristic  function  of  Z  =  X^X^/cr^  a 2  is 


exp 


Y .<*> 


(v-L^+  v^-  2{  v  v2)t2>!  v^i 


it 


2  [  l-(l+p)it  J  [  i+(l-p)it] 


^[l-(l+p)it  J  [  l+(l-p)it  ] 


(H-9) 


Using  properties  of  this  function,  it  is  possible  to  show  that  E 

7 


v^v2+  p  and  the  standard  deviation  is  <T  ■H 


•  V  !P  vi 


V  1  + 1° 


Aroian[  1  Jproved  the  following  statements.  ; 

l)  The  p.d.f.  of  Z  approaches  the  normal  p.d.f.  with  mean  1  and 


variance  <3  as  v,  and  v  - 
z  1  2 


oo  (or  -oo  )  in  any  manner  whatsoever, 


provided  -1  +  £  <^o<l,^.>0. 

2)  The  p.d.f.  of  Z  approaches  the  normal  p.d.f.  with  mean  z  and 


vari'ince  o  ^  if  v, — *  oo  and  v_ 
z  1  2 


■oo 


,  provided  -1<|0  <1-  £, 


£  >  0. 


and  5)  The  p.d.f.  of  Z  approaches  a  normal  p.d.f.  if  v^  remains  constant 


and  v. 


and  v 


-*■  oo  ,  -1  +  £  <  ^  1 1  £  >  0;  or  if  v  remains  constant 

>  -a;  for  -1  <  ^  <T  1  -  £  ,  £  >  0. 


2.  Numerical  Computation  , 

2.1  Integration  of  the  Cumulative  Distribution  Functio n»  The  cumulative 
distribution  function  F(z)  may  be  formulated  directly  by  muling  use  of  the 
fact  that  if  l^(  t  )  is  the  characteristic  function  of  random  variable  Z, 
tnen  the  c.d.f.  of  Z  is  given  by 


oo 


bO 


■W 


cos  t  z 
i  t 


f(-t  )-f(t)  I 


Bin  t  z 


y/(-t)+\f'(t)  dt 


(11-10) 


This  relation  has  been  proved  in[  8  J  and  [  9  ]•  The  advantages  of  this 
formula  lie  in  the  fact  that  a  separate  determin it  ion  of  F(o)  need  not  be 
made  and  a  double  numerical  integration  is  avoided.  Aroian  [2  Juaed  (11-10) 
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to  obtain  numerical  results  when  p  -  0.  In  this  case  F(z)  may  be  expressed 
as 


This  expression  was  numerically  integrated  from  0  to  t^,  t^  to  t^,  .  .  ., 


t0  ti>l»  1  "  l*  2' 


,  where  t^  are  the  zeros  of  sin  |  t 


r  viv2 

L*'  Ut2 

4 


Aroian's  tables  of  this  c.d.f.  include  combinations  of  and  at  inter¬ 
vals  of  0.4,  0  <  v^  <  4,  0  <T  v  <  4.  The  values  of  z  are  given  at  intervals 
of  0.1  for^K^_£  0^}  at  0.2  for  ( ^4  +  <T^ )  to  i  M z+  3  ffj,)  and  for  2“  (5^) 
to  3  (J^)}  at  intervals  of  0.4  for  z+  3  0^)  to  (/^z+  ^  0^)  and  for 

{/U  -  3cr)  to  ( //  -  4  <T  )  and  in  intervals  of  0.8  to  the  extreme  values 
/  z  z  /  z  z 

A- 7  v 

In  theory  the  o.d.f.  of  the  correlated  product  may  be  derived  from 
(II-IO).  However,  the  resulting  expression  is  quite  complicated.  Its 
rather  cumberaome  nature  hinders  the  derivation  of  a  substantial  quantity 
of  numerical  results  using  the  type  of  "intermediate"  computer  dictated  by 
the  scope  of  this  study.  I  sing  ^c^(t)  as  given  by  (II-9),  the  c.d.f.,  F(z), 
of  the  correlated  product  ia 
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2  2  2  2 
wheret  (l-(D  ),  k2-  (l+p  ),  (vJ  +  v2  -  2p  v^). 

In  order  to  obtain  the  zeros  of  the  sin  and  cos  functions,  it  is 
necessary  to  solve  the  fifth-order  polynomial  representing  the  arguments  of 
the  trigonometric  functions.  Numerical  integration  from  zero  to  zero  of 
each  of  these  functions  may  be  accomplished  in  a  number  of  ways  since  all 
derivatives  of  these  trigonometric  functions  are  bounded.  The  number  of 
zeros  of  both  functions  is  greatly  increased  however  in  comparison  with 
(il-ll).  In  addition  a  bound  for  the  tail  area  in  (11-12)  is  difficult  to 
obtain.  Due  to  these  difficulties  and  the  limitations  of  the  available 
computer,  no  further  consideration  was  given  (11-12)  as  a  method  of 
generating  a  large  volume  of  tabular  results. 

2.2.  Integration  to  Obtain  the  Probability  Density  Function!  An 
alternate  approach  t^  obtaining  F(z)  iB  of  course  a  double  numerical 
integration  of  (II-2).  A  rearrangement  of  the  exponents  in  this  equation 
will  allow  z)  to  be  expressed  as 
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4B-/ 


2(1 -p  ) 


z  2  z  /  \ 

-*-(f  vrv2) 

w 


2(l-p  ) 


+  2(p  v.-v^w  ^ 


Vt4VV) 


2(iy) 


dw  , 


(11-15; 


where  K(z)  »  exp  )  - 


2(1^)  L 


r:2+  v22-  2p  z  +  vxv2  J  /2  7T  y~l-/o2  . 


In  turn,  (il-l})  may  be  expressed  as  the  sum  of  the  integrals  of  the  two 


functions  of  the  integrand.  Thu3,  let 
00  00 

i 


K(  z 


j  I1(w)dw  +  /  I2(w)dv  , 


(11-14) 


where  I^(w)  and  I2(w)  are  the  respective  terms  of  ( 1 1—1 3 ) • 

In  order  to  bound  the  tail  areas  of  (II-14)  by  C  ,  L  <  10”^  say,  it 


is  sufficient  to  require  that 


00 


00 


j  Ix(w)dw  +  j  I?(w)dw  <  ♦  £2<C  , 

V  UL„ 


(11-15) 


where  and  axe  the  upper  limits  of  the  numerical  integrations  of 
Ix(w)  and  I2(w),  respectively,  end  and  £2  are  defined  by  (li-16)  and 
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(11-18).  It  Is  possible  to  bound  eaoh  of  these  integrals  by  the  nomal 


probability  integral,  i.e.,  it  is  possible,  for  I^(w)  say,  to  write 


JO 


oo 


j  I^wjdw  -  J  ~ 


•  wife '  * 


2(1-P") 


u. 


< 


=  c 


1 


(11-16) 


The  signs  of  a,(p  v^-v^),  V2~Y1^  combine  to  produce  either 

positive  or  negative  terms  in  the  exponents  of  (11-16).  In  an  error 
analysis,  the  selection  of  an  appropriate  wq  to  insure,  for  all  w  ;>  wq,  that 


(11-17) 


represents  one  method  of  satisfying  (II-16).  Similarly,  the  tail  area 
of  -^(v)  may  be  bounded  by 


dw 
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00 


00 


/ 


I,(v)dw 


<  f 

JTtt  J 

*  M 


2 

w 

2 


dw  -  €, 


(11-18) 


provided  wq  is  chosen  so  that 


exp 


j  i(  z*  '  2  K  VT2)|»  \  f 

I ' 2  V  *2<i-p2)  /  I 


-  < 


•xp 


2  TT 
2 


and 


(H-19) 


for  all  w  >  w  • 
o 

For  all  vq>  y  2  TT  ,  the  sets  of  inequalities  (11-17)  and  (II-19) 

hold  when  the  signs  within  the  exponents  of  these  two  Bets  are  positive. 

In  tnis  oase,  the  upper  Units  of  numerical  integration  of  I^(w)  and 

I2(w)  nay  be  chosen  as  T^«  T?-  4.42,  respectively.  In  this  case  £  £^  and 

froa  an  appropriate  table  each  is  oonputed  to  be  less  than  5  *  10  The 

.5 

total  tail  area  £  is  then  less  than  10  • 

The  other  "extreme"  case  arises  when  the  pair  of  signs  in  either  one 
of  the  sets  of  inequalities,  say  (11-17),  is  negative.  The  first  inequality 
of  ( II— 17 )  then  requires  that 


w  exp 


i  /  *2- 2 1***0  ’rVI* 


IK 


2(.  2s 

w  (1  -p  ) 


)|>  {TW  . 


(11-20) 
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In  most  cases  there  will  be  two  Bets  of  values  of  w  satisfying  this 
inequality.  An  approximation  to  the  least  upper  bound  of  the  upper  range 
of  these  values,  J .  ■ay  be  obtained  by  a  numerical  iteration  of 


2(l-p2)  logj  ~^===.  j  v2- 1 2  z(p  v1-v2)| 


w  +  z  >0. 


(H-21) 


The  set  of  w  satisfying  the  second  inequality  of  (11-17)  is  easily 

12 

jo  v2-vi  |  /(°  •  every  case  the  tail  area  of  (11-16) 


2 

®  .  i- 

1  /  2  -6 
may  be  bounded  by  — ■■ -  /  e  dv  -  £  .  <  5  x  10~  provided  the 

V  2TT  Tx»  4.42 

numerical  integration  is  performed  over  the  interval  (0,  u.  )  where 

L 


a  max 


A  .  2  |  p  v2-v1  |  /p2,  T1=  4.42  (  . 


(11-22) 


The  upper  limit  given  by  (ll-2l)  is  quite  obviously  an  upper  limit 
of  integration  for  I2(w)  by  the  same  argument.  Thus,  the  tail  area  estimate 
of  (11-15)  may  be  restricted  by  £  <  t  ^  £2  <  provided  the  upper 

limits  of  the  numerical  integration  for  both  I^(v)  and  I2(v)  are  determined 
by  (11-22). 


2.3  Methods  of  Numerical  Integration!  Several  formulas  such  as 
Weddel's  formula,  the  trapezoidal  rule,  the  Gregory-Newton  formula,  and  the 
simple  rectangular  formula!"  2  "1,1"  5  "jhave  been  suggested  for  the  numerical 
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integration  of 


I^(w)dw  + 


I2(»)d«. 


(11-25) 


o  O 

The  magnitude  of  the  error  bounds  for  these  and  most  other  numerical 
integration  methods  depends  directly  or  indirectly  on  values  of  a  given 
derivative  of  the  integrand  within  the  interval  of  integration.  The  first 
derivatives  of  I^(w)  and  I^(w)  may  be  written  in  the  form  of  a  rational 
function, 


t  -  Jiw 


Z2(-I)1  +  Iz(p  v1-v2)w-(ln(D2)w2(-l)i+2(p  v2-v1)w5-wi+ 


T. 2v  4 

(1-p  )w 


(11-24) 


where  J^w),  i*»l,  2,  are  the  exponential  functions  of  the  I^(w).  Alter 
nately,  (11-24)  may  be  expressed  aa 


_d_ 

dw 


Ix(w) 


Jt(w) 

(l-fV 


=4 


P(w)  , 


(11-25) 


where  P(w)  represents  the  polynomial  of  (11-24). 

All  derivatives  of  I^(w)  may  be  expressed  in  the  form  (11—25 )  with 
the  order  n  of  the  polynomial  P(v)  increasing  accordingly.  The  terms  of 
P(w)  in  the  second  and  high.r-order  derivatives  are  various  powers  and 
cross-products  of  the  parameters  z,  (p  v2~v^),  (f  ^'^2^  (l-f>  )• 
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The  error  estimate  for  a  given  method  of  numerical  integration  ie 
a  function  of  the  maximum  value  within  the  interval  of  integration  of 
the  derivative  associated  with  that  method.  For  example,  the  error  bound 
of  the  trapezoidal  rule  is  a  function  of  the  second  derivative,  Weddel's 
formula  involves  the  sixth  derivative,  etc.  The  actual  error,  of  course, 
may  be  much  smaller  than  indicated  by  the  error  bound. 

In  order  to  calculate  the  maximum  values  of  the  given  n®  derivative 
within  the  interval  of  integration,  it  is  necessary  to  derive  and  solve  for 
the  roots  of  the  polynomial  of  the  (n+l)8*  derivative.  This  task  becomes 
increasingly  difficult  to  do  analytically  as  the  order  of  P(w)  ir  -reaaes. 

A  numerical  iteration  method  must  be  used  to  solve  for  the  roots  of  P(v)  in 
the  higher-order  derivatives. 

The  values  of  each  derivative  are  functions  of  the  values  or  the 
parameters  and  their  signs.  Thus,  for  appropriate  sets  of  parameter  values, 
the  derivatives  are  large  in  the  neighborhood  of  w  =  0,  or  more  generally 
in  the  interval  (0,l).  As  yet,  no  method  of  characterizing  the  derivatives 
within  the  interval  of  integration  as  functions  of  the  parameter  values  has 
proved  satisfactory.  Because  of  these  difficulties,  other  numerical 
methods  which  are  not  functionally  dependent  of  the  derivatives  are  believed 
to  be  more  expedient  for  this  problem. 

The  c.d.f.'s  of  z  appearing  in  Appendix  B  were  obtained  using  a  double 
numerical  integration  of  (11-23)  by  the  simple  upper  and  lower  sum  rectangular 
formula.  In  order  to  obtain  ^ (z),  the  real  roots,  R^,  R^  .  •  . ,  of  P(w) 


2fl 


for  the  first  derivative  of  both  I^(w)  and  (w) #  ( 11-24 ) ,  were  estimated 
by  numerical  methods  to  within  five  significant  digits.  For  eaoh  integrand 
Ii(w),  an  upper  and  lower  sum,  t\::  and  L^,  were  computed  in  the  intervals 
(0  to  ,  (Rj  to  U^)  with  a  normal  increment  Av  .  0.01.  The  increment 

was  reduced  as  necessary  to  insure  that  Ug-Lg<  10  in  all  intervals.  In 
many  cases,  the  functions  I^(w)  are  quite  steep  in  the  interval  (0,  R^)  with 
R^<  <  1.  In  these  oases,  very  small  A  w's  were  required  to  obtain  a 
satisfactory  estimate  in  this  interval. 

As  compared  with  other  numerical  methods,  the  rectangular  formula 

provides  greater  accuracy  but  generally  requires  a  much  larger  number  of 

computer  calculations.  The  actual  computer  time  required  to  obtain  F(z) 

is  dependent  upon  the  shapes  of  I^w),  I2(w),  and  <^?(z).  The  "average" 

computer  time  required  to  approximate  the  double  integration  within  the 

desired  accuracy  was  approximately  5W  minutes^.  The  computer  time  required 

for  this  integration  program  can  be  reduoed  to  approximately  1.25  minutes 

4 

using  a  high-speed  computer  such  as  an  IBM  7094  ,  Considering  both  com¬ 
puter  costs  and  the  time  required  to  generate  the  desired  volume  of 
tabular  results,  the  use  of  an  IBM  7094  or  its  equivalent  is  recommended 
for  future  work.  The  cost  of  using  a  computer  on  this  scale  prohibited  its 
use  in  the  investigation  which  was  intended  and  funded  as  a  "preliminary 
study". 


3 

4 


Based  on  the  use  of  on  IBM  1620  computer. 

As  estimated  by  members  of  the  Martin  Company  Data  Systems  Division 
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appendix  a 


We  give  two  results  here  that  were  obtained  in  a  "side"  investi¬ 
gation  involving  the  quotients  of  random  variables.  One  is  probably 
not  new  but  an  easy  proof  is  given. 

n 


Theorem  Let  f(z)  -  ^  ^  b^  *n  *  be  a  oomplex  polynomial,  bQ  ^  0. 


1 — O 


Then  the  zeros  of  f(z)  are  in  disk  z  <  1  +  where  p  -  sup 


i  a  1 1  •  •  •  i  n  o 


Proof!  Write  f(z)  =  b  g(z)  where 


n 


«(«) 


-  z  + 


Y  ai  i“'1 


i-i 


and 


a 

l  b 

o 


Then 

|  g(z)| 


v  i-1  ' 


n 


>  * 


Y ai ini  -  I1"! -X^i°iiNn'i- 


i.l 


i-1 


Let  p  -  sup  a^|  and  suppose  Jzj>  1  ♦  p,then  from  the  above 


M 

g(z)|>|z|n-p  ^2  [zl"-1  -  |z|p-  p 


i.l 


In  T 

z  -  1 


*  -  1 


z 

n  ( 

1  3 

*  e 

z 

-  1 

But  |  s  {  -  1  “f>>  1  +  p  -  1  -  p  -  0  and  hence  for  |  s  |  >  1  ♦  pt  j  ^(z)  |  o 
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Consequently  |f(z)|  -  j  bQ  |  |g(z)|>0  for  thiu  range  of  z  and  the  result 
follows. 

The  question  arose  in  a  study  of  mappings  between  various  topo¬ 
logical  spaces,  wnen  the  arbitrary  union  of  closed  sets  is  closed.  As  a 
partial  answer  we  give  the  following.  Notationally  for  the  space  Y,  let 

2Y  ECy(  E  is  closed  and  non-empty j  .  A  mapping  f,  from  a  space  X  into 
Y 

2  is  said  to  be  upper  semi-continuous  in  case  xq£  X,  U  open  in  Y  and 
f(xQ) d  U  implies  that  there  is  an  open  set  V  in  X,  xq£  V,  such  that  x  £  V 
implies  that  f(x) C  l  the  above  to  hold,  of  course,  for  each  xq£  X. 

Theorem!  Let  X  be  a  compact,  Hausdorff  space,  Y  regular,  and  let  f  be 


an  upper  aemi-continuous  function  from  X — *•  2  . 
is  closed  in  Y. 


Then 


U 

x£  X 


f(x) 


Proof i  Let  B  -  y'  f(x)  and  assume  y  £  B,  y  £  B  where  B  denotes  the  closure 

of  B  in  Y.  Clearly  then  y  f(x)  for  any  x.  This  implies  that  for  each 

x  there  are  two  disjoint  open  sets  in  Y,  and  ®yf(x)  8UC^  ^hat 

f(x)c  Wf^  and  y  £  ^(x)*  Qyf(x)^  Wf(x)  "  ^  3ince  f  is  upper  semi" 

continuous,  for  each  x  €  X  there  is  an  open  in  X,  x  £  V^,  such  that 

f(x) <z  ^or  jach  x  £  V^.  Hence  jv^j  is  a  cover  for  X,  a  compact  space. 

n 

Consequently  a  finite  number  will  do*  say  X  =  V  •  Then 

i-1  xi 

n  n 


H  2  -  U  W.  ,  «  3  U  f(x)  -  B.  Now  0  3  f\  0  /  \  is  an  open  set 

i_l  x£  X  »  i^l  ** 

containing  y,  and  0  D  Z  =»  0  so  that  0  f)  B  *  0.  But  this  is  a  contradiction 

y  y 

since  y  was  assumed  to  be  in  B  and  the  theorem  is  proved. 
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APPENDIX  B 

TABLES  OF  THE  PRODUCT  OF 
TWO  NORMALLY  DISTRIBUTED 
RANDOM  VARIABLES 
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1#  The  c.d.f.'s  of  the  random  variable  Z  -  for  various 

parameter  values  were  obtained  by  a  double  numerical  integration  of  (II-2J) 
In  this  preliminary  study,  only  positive  parameters,  ^  ,  v^,  v^,  were 

considered. 


2.  Tail  Area  Bounds.  Denoting  the  p.d.f.  of  Z  in  the  correlated  case 

as  f (z  ,  ^>0),  it  is  easily  shown  that  f(z#  ^  >  0)<  <p  (z  ,  ^  0)  for 

z  z 

z  <0.  Thus,  r  f  (z,  ^  >  0)  dz  <  f*  (*.  S-0)  dz  «  A  v  for  Zq  <  0. 


-  CD 


00 


some  z^>0,  it  follows  that 


The  value  z^  may  be  chosen  so  that  A^  is  arbitrarily  small  .  In  addition, 
it  may  be  shown  that  f(z,  ^>0)  <  (p^max  |  (Xj/  ®i>?  (*2/  <T2)2|)  for  all  z  >  z 

Here  the  symbol  <j)  ^max  |  (xl/  (*2/  <T>)  1)  denotes  the  p.d.f.  of  the 

square  of  the  largest  of  the  random  variables  x^/  <T^  and  x^/  <T,.  This 

2 

random  variable  follows  the  x  p.d.f.  with  one  degree  of  freedom.  Thus  for 

f(z,  ^>0)dz  max  |  (x^/ 

(x2/  *,)  1)  i(  •  )  “  A^«  The  values  zq  and  z^  may  be  chosen  so  that  the  sum, 
0,  of  the  probabilities  A^  and  as  determined  by  their  respective  c.d.f.'s 

is  arbitrarily  small.  The  integral  value  I  of  f(z  ,  ^  )  for  positive  ^ 
may  be  estimated  for  the  neighborhood  (  <  0<A^)  containing  the  point 

of  discontinuity,  z=0,  from  the  relation 

I  =*  1  -  |  f(z  ,  ^>0)dz  +  f  1  f(z,  q>0)dz  V  e). 

'  Zo  *7  j 


The  estimate  I.  may  be  accurately  aj-proximaU  d  simply  by  requiring  that  0 
be  made  small. 


* 

This  function  has  been  tabulated  by  L.  A. 


Aroian  [  2  ]  . 


>0. 
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3.  Checks.  Tables  B-I  to  B-mmay  be  oompared  with  Aroian's  results  [?  J  . 

-5 

Those  data  points  noted  with  an  asterisk  vary  by  10  with  his  results. 
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B  -  I  . 

Parameter  Values; 

?  .  0.,  Vj  .  0.,  V2  .  0.4 


z 

f(*) 

r(z) 

z 

f(z) 

F(z) 

-7.60 

.00014 

.00011 

.10 

82755 

.60294 

-6.80 

00028 

.00025 

.20 

.55167 

.66551 

-6.00 

.00064 

.00057 

•  30 

.43129 

.71327 

-5.20 

.00084 

.00128 

.40 

•  35191 

■75177 

-4.40 

.00287 

.00290 

•  50 

.29379 

•  78369 

-4.00 

.00468 

.00439 

.60 

.81053 

-3.60 

.00715 

.00665 

•  70 

.21301 

.83344 

-3.20 

. 01027 

.01011 

.80 

•19370 

.85313 

-3.00 

.01332 

.01249 

•  90 

.15935 

. 87018 

-2.80 

.01658 

.01544 

1.00 

.13489 

.88500 

-2.60 

.02067 

.01912 

1.20 

.10770 

.90931 

-2.40 

.02854 

. 03371 

1.40 

.08356 

.92809 

-2.20 

.03240 

.02946 

1.60 

.06560 

•94273 

-2.00 

.04076 

.03667 

1.80 

.05171 

.95424 

-1.80 

.05171 

.04576 

2.00 

. 04076 

.96333 

-1.60 

.06560 

.05727 

2.20 

.03240 

.97654 

-1.40 

.08356 

.07191 

2.40 

.02584 

.97629 

-1.20 

.10770 

.09069 

2.60 

.02067 

.98088 

-1.00 

.13489 

.11500 

2.80 

.OI658 

.98456 

-  .90 

.15925 

.12982 

3-00 

.01332 

. 98751 

-  .80 

.18370 

.  14687 

3-20 

.01027 

.98989 

-  .70 

.21301 

.16656 

3.60 

.00715 

•99335 

-  .60 

.24888* 

.18948 

4.00 

.00468 

.99561 

-  .50 

.29379 

.21634 

4.40 

.00287 

•  99710 

-  .40 

. 35191 

.24823 

5.20 

.00145 

.99872 

-  .30 

.43129 

.28672 

6.00 

.00064 

.99943 

-  .20 

.55167 

•33449 

6.80 

.00028 

.99975 

-  .10 

.82755 

.39706 

7.60 

.00014 

.99989 

.00 

00 

.50000 

L 
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B  -  II  . 

Parameter  Values} 

^  -  0.  ,  V  -  0.4,  V2  -  0.4 


z 

f(z) 

F(z) 

Z 

f(0 

F(z) 

.00005 

.00004 

.20 

.00011 

.00010 

•30 

.65678  | 

.00026 

.00022 

.40 

.69828 

.00063 

.00053 

.50 

.  31109 

. 73173 

.00129 

.00124 

.60 

.26905 

.76050 

.00216 

.00192 

.70 

.23494 

.78554 

-4.00 

.00337 

.00297 

.80 

.20661 

.00749 

-3.60 

.00530 

.00462 

•  90 

.18268 

-3.20 

.00777 

.00721 

1.00 

. 16221 

.84403 

-3.00 

.01024 

.00903 

1.10 

.1^453 

.85930 

-2.80 

.01292 

.01131 

1.20 

•  12915 

.87293 

-2.60 

.01634 

.01420 

1.30 

.11569 

.88513 

-2.40 

.02070 

.01785 

1.40 

.10135 

.89607 

-2.20 

.02633 

.02248 

1.60 

.08453 

.91474 

-2.00 

•03359 

. 02838 

1.00 

.06874 

.92989 

-1.00 

.04303 

.03592 

2.00 

.05615 

.94224 

-1.60 

.05543 

.04556 

2.20 

.04600 

.95234 

-1.40 

.07187 

.05809 

2.40 

•03775 

.96063 

-1.20 

.09398 

.07434 

2.60 

.03106 

.96744 

-1.00 

.11913 

.09568 

2.80 

.02560 

.97305 

-0.90 

.14244 

.10884 

3-00 

. 02118 

.97768 

-0.80 

. 16476 

. 12407 

3-20 

.01746 

.98151 

-0.70 

.19295 

.14 180 

3-40 

.01444 

.98467 

-0.60 

. 22646 

. 16256 

3.60 

.01148 

.98728 

-0.50 

.26920 

.18709 

4.00 

.00834 

.99124 

-0.40 

•  32465 

. 21640 

4.40 

.00574 

.99396 

-0.  30 

.40048 

.25204 

4.80 

.00370 

.99583 

-0.20 

•51533 

.29650 

5.60 

.00201 

.99801 

-0.10 

.  77595 

•35509 

6. 4o 

.00095 

•99905 

0.00 

00 

.45169 

7.20 

.00045 

•99955 

.10 

. 79584 

.54958 

i  8.00 

.00028 

•99978 

j 

39 


B—  III 


Parameter  Values; 


•11.00 

10.00 

9.U0 

<J.6o 

7-80 

7.00 

6,20 

5.80 

5.40 
5.00 

4.60 

4.20 
4.00 
3.80 

3.60 
3-40 

3.20 
3-00 
2.00 
2.60 

2.40 

2.20 
2.00 
1.80 
1.60 

1.40 
1.20 
1.10 
1.00 

•  90 
.80 

•  70 
.60 

•  50 
.40 

•  30 
.20 
.10 
.00 
.10 
.20 


*  =  0.,  Vx  =  0.4,  V2  =  1.2 


f(r) 


.00001 

.00002 

.00004 

.00009 


.00827 

.01059 

.01279 

.01548 

.01874 

.02272 


.07666 
.09096 
.10414 
.11659 
.13589 
. 14748 
.16702 
.19030 
.21865 
•25425. 
.30128 
.36983 
.51705 
00 

.55191 

.42497 


.00001 

.00002 

.00004 

.00009 

.00021 

.00042 

.00084 


.00369 

•00535 

.00644 

.00775 

.00934 

.01126 

.01358 

.OI638 
.01977 
.02388 
.02886 
.03492 
.04229 
.05129 
.06229 
.07583 
.09256 
. 10239 
.11338 
.12570 
•13955 
.15519 
.17295 
.19325 

.21669 

.24411 

.27694 

.31807 

.38035 

.44400 

.49073 


•  30 
.40 

•  50 

.60 

•  70 
.80 

•  90 
1.00 
1.10 
1.20 
1.30 

1.40 
1.50 
1.60 
1.70 
1.00 
1.90 
2.00 
2.10 
2.20 

2.40 
2.60 
2.80 
3.00 
3-20 

3.40 
3.60 

3.80 
4.00 

4.20 

4.40 

4.60 

4.80 
5.00 

5.20 

5.40 

6.20 
7.00 

7.80 

8.60 

9.40 


.29967 
.27002 
. 25287 
.23380 
.21672 
.20116 
.18775 
. 17418 
.16223 
.15116 
.14089 
.13137 
.11999 
.11170 
.10649 
.09930 
.09258 
.08489 
.07517 

.06527 

.05666 

: 84256 
.03685 
.03244 
.02758 
.02383 
.02058 
.01776 
.01531 

-m 

.00974 

.00881 

.00594 

.00263 

.00141 

.00075 

.00040 


F(z) 


.53983 
•56451 
•59584 
.634  44 
.65074 
.67502 

.69750 
.71846 
•  73775 
.75579 
.77259 
.78824 
.80282 
. 81642 
.82909 
.84091 
.85195 
.86221 
.87180 
.88073 
.89682 
.91000 
.92293 
.93345 
•94257 
. 95048 
.95732 
.96323 
.96835 
.97277 
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B-  IV. 

Paraaeter  Values; 


*  =  o.t  vl=  0.5,  v2=  0.5 


r(z) 


f(z) 


-7-20 

•6.40 


.00001 

.00002 

.00004 

.00010 

.00023 

.00054 

.00128 

.00311 

.00408 

.00769 

.01221 

.01954 

.03164 

.05123 


.00022 
.00053 
.00126 
.00302 
..  00469 
.00717 
.OLIOS 
.01743 

.02765 

. 04424 
07210 
.  12087 
.21286 
.28465 
.34084 
.42666 


.10 

.20 

.40 

.80 

1.20 

1.60 

2.00 

2.40 
2.80 
3.20 

3.60 
4.00 
4.80 

5.60 

6.40 


.68106 
.  51947 

.36322 

. 21768 
.14230 
.09621 
.08589 
.06614 
.04589 
.02241 
•01571 
.01103 
.00545 
.00269 
.00133 


r(z) 


. 47012 
. 53312 

.61139 

.72757 

•79956 

.84726 

.88368 

•91409 

.93650 

•95116 

•95778 

.96313 

.96815 

.97206 

•97467 


8S  setts'?  S3  @8S»?$S'*,S,'#>S<j8£8"#. 


B-  V. 

Parameter  Values; 

=  0.3,  V.=  0.15,  V  =  0.4 


f(*) 


P(*) 


-6.84 

-6.12 

-5-40 

-4.68 


.00001 

.00003 

.00008 


12 
.01504 
.02017 
.02729 
.03690 
.05012 
06887 

.09195 
.11209 
.13350 
.15904 
.19203 
.23503 
.29068 
.36784 
.48580 
. 75244 


.00002 

.00006 

.00010 

.00014 

.00041 


•00537 

.00673 

.00773 

.00950 

.01237 

.01523 

.01953 

.02525 

.03309 

.04382 

.55739 

,06149 

.07762 

.09082 

.10637 

.12593 

.14957 

.17924 

.21765 

•27337 

.36359 


•  09 
.18 
.27 
.36 
.45 
.54 

.63 

.72 

.81 

•  90 
1.08 
1.26 
1.44 
1.62 
1.80 
1.98 
2.16 
2.34 
2.52 
2.70 
2.88 
3-24 
3.60 
3.96 
4.68 
5.40 
6.12 
6.84 


.80128 
. 55214 
.44570 
•  37549 

.32367 
.28311 
.25018 
.22277 
•19953 
.17439 
»345 
. 12278 
. 10276 
.08636 
.07257 
.06149 
.05229 
.04460 

.03813 

.03267 

.02684 

.02125 

.01500 

.01101 

.00723 

.00410 

.00234 

.00154 


.45604 

. 51638 
.56189 

.59884 
.62030 
.65760 
.68179 
•70237 
.72189 
. 73872 
75761 
.79271 
.81239 
.82975 
.85317 
.85597 

. 86603 

.87371 

.88228 

nAn// 


.89401 

.89761 

.90933 

.91415 

.91744 

.92039 

.92286 

.92831 
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B-  VI 


Paraaeter  Valuea; 


^  =  0.3,  V1=  0.55,  V2=  0.55 


z 

r(z) 

r(0 

Z 

r(*> 

Hz) 

-7.18 

.00000 

.00000 

.09 

.72204 

•37697 

-6.46 

.00000 

.00000 

.18 

.51147 

.44612 

-5.71* 

.00003 

.00000 

•  27 

.42280 

.48066 

-5.02 

.00009 

.00001 

.36 

.  36426 

.52420 

-4.30 

.00024 

.00002 

•  45 

.32079 

.56513 

-3.94 

.00046 

.00017 

.54 

.28645 

•58245 

-3-59 

.00082 

.00039 

.63 

.25827 

.60697 

-3.23 

.00147 

.00079 

•  72 

.23451 

.62914 

-2.87 

.00245 

.00134 

.80 

.21409 

.65110 

-2.69 

.00344 

.00247  1 

.89 

.19629 

.66780 

-2.  C1 

.00463 

.00302 

.98 

•  18058 

.67032 

-2.33 

.00624 

00391 

1.07 

.16661 

.68595 

-2.15 

.00844 

.00490 

1.16 

. 15410 

.70038 

-1.97 

.01144 

.OO669 

1.25 

•  13939 

.71359 

-1.79 

.01556 

.00811 

1.43 

.12395 

.72865 

-1.6l 

.02125 

.01243 

1.6l 

.10746 

.74938 

-1.43 

.02919 

.OI696 

1.79 

.09358 

.77621 

-1.25 

.04034 

.02323 

1.97 

.08172 

.79198 

-1.07 

.05624 

.03190 

2.15 

.07150 

.80578 

-  .89 

.07601 

.04382 

2.33 

.06272 

.81786 

-  .80 

.09384 

.05146 

2.51 

.05511 

.82837 

-  .72 

.11207 

.06073 

2.69 

.04860 

.83779 

-  .63 

.13552 

•07187 

2.87 

.04271 

.84601 

-  .54 

.16422 

.08436 

3.24 

03093 

.85926 

-  .45 

.20157 

. 10181 

3-59 

.02062 

.86854 

-  .36 

.25100 

.13219 

4.30 

.01511 

.88016 

-  .27 

.31969 

. 14787 

5.02 

.01060 

.89066 

-  .18 

.42475 

.18137 

5.74 

.00652 

.89682 

-  .09 

.66035 

.23020 

6.46 

.00402 

.90061 

.00 

.297^7 

7.18 

.00320 

.90322 
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B-  VII. 

Parameter  Values; 


^  =  0.3,  V1=  0.7,  V2=  0.7 


z 

f(z) 

- - 

F(*) 

z 

f(z) 

F(z) 

-8.6o 

.00000 

.00000 

-  .18 

•  37*09 

.15722 

-7.90 

.00000 

.00000 

-  .09 

•53018 

.19792 

-7.20 

.00000 

.00001 

.00 

00 

. 26839 

-6.45 

.00000 

.00002 

.09 

. 59228 

•  3**l65 

-5.75 

.00002 

.00003 

.18 

. 46645 

.38929 

-5.00 

.00007 

.00004 

.36 

•  3^770 

.46256 

-4.30 

.00023 

.00006 

.72 

. 23683 

.56770 

-3.60 

.00072 

.00009 

1.08 

•  17597 

.64208 

-3.25 

.00129 

.00437 

1.44 

.13521 

.69809 

-2.88 

.00232 

.00727 

1.80 

.13720 

. 74721 

-2.52 

.00419 

.00844 

2.16 

.12008 

•  793*0 

-2.16 

.00763 

.01057 

2.52 

.09470 

.83107 

-1.80 

.01405 

.01447 

2.88 

.05256 

.85847 

-1.44 

.02586 

.02165 

3.25 

.04187 

.87560 

-1.08 

.05052 

.035^0 

3.60 

.033*0 

.88908 

-  .72 

. 10151 

.06277 

4.30 

.02135 

.89105 

-  .36 

.22545 

.12163 

5.00 

5-75 

_ — J 

.01363 

.00871 

.89214 

.90100 
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B-  VIII. 

Parameter  Values; 

^  *  0.45,  V1=  0.2,  V2=  0.45 


2 

r(0 

?(z) 

z 

f(z) 

r(z) 

-6.10 

.00000 

1 

.08 

.76179 

.44780 

-5.46 

.00001 

.00001 

.16 

.53068 

.49950 

-4.82 

.00004 

.00002 

.24 

.43354 

.53807 

-4.18 

.00023 

.00010 

.32 

.36966 

.57175 

-3.52 

.00036 

.00029 

.40 

.32249 

.59788 

-3-20 

.00070 

.00045 

.48 

. 28548 

.62220 

-2.88 

.00129 

.00078 

•  56 

.25533 

. 64  383 

-2.56 

.00221 

.00134 

.64 

.23009 

.66325 

-2.40 

.00313 

.00185 

•  72 

.20858 

.68079 

-2.24 

.00426 

.00236 

.80 

.18450 

.69652 

-2.08 

.00580 

.00316 

.96 

.16087 

-1.92 

.00874 

.00432 

1.12 

.13628 

.71315 

-I.76 

.01084 

.00589 

1.28 

.11684 

.73340 

-1.60 

.01489 

.00795 

1.44 

.10057 

.75080 

-1.44 

.02063 

.01079 

1.60 

.08656 

.76577 

-1.28 

.02858 

.01473 

1.76 

.07514 

.81347 

-1.12 

.03975 

.02019 

1.92 

.06544 

.82472 

-  .96 

.05595 

.02785 

2.08 

.05717 

.83453 

-  .80 

.07652 

.03845 

2.24 

.05007 

.84310 

-  .72 

.09441 

.04529 

2.40 

.04394 

.85063 

-  .64 

.11380 

.05362 

2.56 

.03698 

.85710 

-  .56 

.13790 

.06368 

2.88 

.C3071 

.85795 

-  .48 

. 16837 

.07593 

3.20 

.02397 

.85806 

-  .40 

.20769 

.09098 

3.52 

.01752 

.85972 

-  .32 

.25997 

.10968 

4.18 

.01266 

.85997 

-  .24 

.33295 

.13349 

4.82 

.00791 

.86003 

-  4 16 

.44504 

.16452 

5.46 

.00498 

.86117 

-  .08 

.69762 

. 21023 

6.10 

.00361 

.86203 

.00 

CO 

.31104 
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B-  IX 


Parameter  Values; 


< 

£  =  0.45,  V1=  0.8,  V2= 

0.8 

f(e) 

r(«) 

z 

f(z) 

F(z) 

-7-75  .  00000  .  00000  -  .08  .1*5891  .15255 

-7.10  .00000  .00000  .00  ©  .20720 

-6. 45  .  00000  .  00000  .  08  .  52512  .  26550 

-5-do  .00000  .00000  .16  .1*1855  .30224 

-5.15  .00001  .00001  .32  .31957  .36129 

-4.50  .  00003  .  00002  .  64  .22838  .41*751 

-3.85  .00011  .00006  .96  .17803  .51239 

-3.22  .  00039  .  00022  1.28  .14353  56398 

-2.90  .00073  .00041  1.60  .15279  -61139 

-2.58  .00138  .00073  1-92  .14030  .65329 

-2.25  .  00262  .  00128  2.25  .11609  69931 

-1.92  .00501  .00259  2.58  .06760  .72873 

-1.60  .00968  .00492  2.90  .05650  .74857 

-1.28  .01369  .00946  3-22  .04733  .76518 

-  .96  .03832  .01860  3.85  .03327  .79137 

-  .64  .  08078  .  03766  4.50  .  02339  .  80910 

-  .32  .18825  .  08071  5.15  .  01644  .  82184 

-  .16  .32020  .12139 


B-  X 


Parameter  Values; 


$  =  0.5,  v1=o.f 


V  0. 


f(t) 


r(z) 


-5.250 

-4.650 

-4.050 

-3.450 

-2.850 

-2.550 

-2.250 

-2.100 

-1.950 

-1.800 

.1.650 

1.500 

1.350 

1.200 

1.050 

.900 

.750 

.675 

.600 

.525 

.450 

.375 

.300 

.225 

.150 

.075 

.000 


.00000 

.00001 

.00005 


.00001 

.00003 


i 


.00205 

.00297 

.00408 

.00528 

.00827 

•01155 

.01635 

.02328 

.03335 

.04818 

•07039 

.08555 

.10446 

.12831 

.15878 

.19846 

.25154 

.32713 

.43717 

.63642 


.00083 

.00105 

.00173 

.00210 

.00259 

.00332 

.00417 

.OO569 

.00764 
.01049 
.01401 
.02924 
.03217 
.04554 
.05816 
.07253 
•09039 
.11238 
. 14183 
.17004 
.18215 
.23574 


z 

f(z) 

rM 

.075 

.70335 

.49978 

.150 

•53396 

.54327 

.225 

.44158 

.58262 

•  300 

•37526 

.61370 

.375 

.32720 

.63959 

.450 

.28931 

.66271 

.525 

.25838 

.68325 

.600 

.23250 

.70166 

.675 

.21042 

.71837 

.750 

.19135 

.73347 

.900 

•15997 

.76391 

1.050 

.13526 

.78235 

1.200 

.11531 

.80114 

1.350 

.09895 

.81683 

1.500 

.08535 

.83021 

1.650 

.07469 

.84211 

1.800 

.06426 

.85208 

1.950 

.05499 

.86203 

2.100 

.04899 

.86983 

2.250 

.04134 

.87660 

2.550 

.03432 

.88372 

2.850 

.02477 

.89684 

3.450 

.01793 

.93846 

4.050 

.01108 

95717 

4.650 

.00692 

.96259 

5.250 

.00490 

.96653 
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B-  XI 


Parameter  Values; 


^  =  0.5,  V  =  0.2,  V2=  0.45 


z 

f(z) 

F(z) 

z 

f(0 

F(z) 

-5.700 

.00000 

.00000 

.075 

.73360 

. 39832 

-5.100 

.00000 

.00000 

.150 

.51411 

. 4  3256 

-4.500 

.00001 

.00001 

.225 

.42254 

. 50764 

-3.900 

.00002 

.00002 

.  300 

. 36244 

•53708 

-3.300 

.00026 

.00004 

.375 

. 31810 

. 56260 

-3.000 

.00053 

.00026 

.450 

. 28329 

•57515 

-2.700 

.00099 

.00039 

.525 

.25489 

.60533 

-2.400 

.00174 

.00080 

.600 

. 23108 

.62357 

-2.250 

.00250 

.00122 

.675 

. 21074 

.64013 

-2.100 

.00344 

.00156 

.750 

.18753 

.65506 

-1.950 

.00475 

.00218 

.900 

.16448 

. 67626 

-1.800 

.00724 

.00308 

1.050 

.14189 

.70031 

-1.650 

.00909 

.00403 

1.200 

.12311 

.72019 

-1.500 

.00593 

1-350 

. 10725 

. 73646 

-1.350 

.01772 

.OO836 

1.500 

.09342 

.75151 

-1.200 

.02485 

.01142 

1  . 65O 

.08207 

•75467 

-1.050 

.03498 

.01587 

1.800 

. 07234 

. 77626 

-  .900 

.04984 

.02284 

1.950 

.06379 

-  .750 

.06899 

.03139 

2.100 

.05670 

.80458 

-  .675 

. 08562 

.03321 

2.250 

.05036 

-  .600 

.^0383 

.03679 

2.400 

.04289 

.81583 

-  .525 

. 12657 

.04543 

2.700 

.03648 

.83106 

-  .450 

.155^7 

.05601 

3.000 

.02917 

-  .375 

.19293 

.06809 

3.300 

.02184 

. 84665 

-  .300 

.24295 

.0b};2 

3.900 

.0 1656 

.85117 

.  .225 

.  31302 

. 10627 

4.500 

.01085 

-  .150 

.42092 

.13369 

5-100 

.00717 

. 87279 

-  .075 

.66379 

.17344 

5.700 

.00546 

•  87859 

.OX) 

00 

.26058 

B-  XII. 


Parameter  Values; 


-6.000 
-5-400 
-4.800 
-4.200 
-3.600 
-3.300 
-3.000 
-2.700 
-2.400 
-2.250 
-2.100 
-1.950 
-1.000 
-I.65O 
-1.500 
-1.350 
-1.200 
-1.050 
-  .900 

-  .750 

-  .675 

-  .600 

-  .525 

-  .450 

-  -375 

-  .300 
.225 
.150 
.075 
.000 


.00002 

.00008 


.00064 
.00115 
.00168 
.00235 
.00328 
.00460 
.00646 
.00811 
.01291 
.OI838 
.02633 
.03806 
.05332 
.06702 
.08150 
.10034 
.12380 
.15471 
.19615 
.25438 
•  34411 
.54470 
00 


0.5,  V1=  0.7,  v?=  0.7 


r(z) 


.00000 

.00000 

.00001 

.00001 

.000C7 


.00577 

.00610 

.00677 

.00762 

.00879 

.01043 

.01279 

.01594 

.02312 

.04006 

.05457 

.06015 

.06697 

.07537 

.08582 

.09895 

.11586 

.13829 

.17162 

.21884 


•  075 
.150 

•  225 
300 

•  375 
.450 
.525 
.600 
.675 

•  750 
.825 
.900 

•  975 
1.050 
1.200 
1.350 
1.500 
1.650 
1.800 
1.950 
2.100 
2.250 

2.400 
2.550 
2.700 
3.000 
3.300 
3.600 
4.200 
4.800 

5.400 
6.000 


f(t) 


.32876 
.29478 
.26801 
.24603 
.22746 
.21143 
•19737 
.18  487 
.17366 
.16354 
.15062 
.13884 
. 12478 
.11265 
.10198 
.09250 
.08412 
.06662 
.07005 
.06382 
.05833 
.04877 
.04552 
.03823 
.03012 
.02440 
. 01734 
.01234 
.01136 


r(*) 


.29582 

•  33567 

•  35643 

.38280 

.40161 

.43729 

.45657 

.47369 

.49007 

.50611 

.52061 

.53389 

. 54654 

.55832 

.57619 

.60040 

.61760 

.63370 

.64829 

.66154 

.67112 

.68459 

.69463 

•70379 
.71183 
. 72146 
.72998 
•73513 
•75149 
.76401 
.77290 
.78003 
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B- XIII 


Parameter  Valueo; 

^  =  0.5,  v1=  0.85,  v2=  0.85 


z 

r(t) 

F(z) 

z 

f(z) 

F(z) 

1 

.00000 

.00000 

-  .150 

.28920 

.09971 

I 

.00000 

.00000 

-  .075 

.41698 

.15167 

.00000 

.00000 

00 

•  17293 

-5.400 

.00000 

.00000 

.075 

.48290 

.22214 

-4.800 

.00000 

.00000 

.150 

. 38721 

.25477 

-4.200 

BWrrH 

.00002 

.  300 

. 29921 

. 30625 

-3.600 

.00003 

.600 

.21902 

. 38216 

-3-000 

UrrrraJ  ■' 

.00014 

.900 

. 17488 

.44297 

-2.700 

.00026 

1.200 

.14441 

.49086 

-2.400 

.00104 

.00052 

1.500 

.157U7 

•53617 

-2.100 

.00203 

.00099 

1.800 

. 14811 

. 58198 

-1.800 

.00397 

.00198 

2.100 

•12553 

. 62314 

-1.500 

.00785 

.00362 

2.400 

.07488 

.65319 

-1.200 

.01552 

.00714 

2,700 

.06410 

. 67404 

-  .900 

.03259 

.01434 

3.000 

.05499 

. 69191 

-  .600 

.07038 

.02980 

3-600 

. 04056 

.72056 

-  .300 

.16799 

.06556 

4.200 

.02992 

.74171 

4.800 

.02207 

.75745 

-  . 
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